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(71) We, Rolls-Royce Limited a British Company of Moor Lane 3 Derby, 
do hereby declare the invention, for which we pray that a patent may be granted to 
us, and the method by which it is to be performed, to be particularly described in and 
by the following statement: . . t e 

5 This invention relates to an optimizing control system. . -> 

According to this invention there is provided a system for varying an input or a 
manufacturing or power plant automatically to drive an output of the plant to a 
maximum or minimum value, comprising a digital computer adapted to act on said 
input and respond to said output and programmed to perform iteratively the opera- 
tions of converting a pseudo-random binary sequence into a perturbation signal where* W 
in one of the two types of digit of the sequence produces a ramp signal component 
successive ones of which are of opposite slope and the other type of digit produces 
signal components of uniform level connecting the successive ramp components, appiy- 
in« the perturbation signal to the plant input, reading the consequential perturbations 
oAhe plant output^ correlating the perturbation signal and the consequential perturba- 
tions to determine the slope of the innut/outnnr «»io*:^-w- 
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25 THE PATENT OFPICE 

* ■* v 

r A#0 tt uuiuct- xh 3 a turbine 16 adapted to drive the compressor, a 

jet tube 18 *and a variable area exhaust nozzle 20 supported by the jet tube. The 
burner is connected to a fuel supply controlled by a valve 22 which is operable by 
a signal 1 The nozzle 20 is connected to an actuator 24 for varying the area of the 30 
nozzle. The actuator 24 is operated by a signal 2. The signals l a 2 are hereinafter also 
referred to as "parameters". 

Variation of the parameters 1, 2 causes changes in the output X of the engine as 
measured by a thrust sensor 26. The variation of X is shown by the characteristic 
35 in Fig 2 where it is defined by the area 12X. The maximum or optimum thrust 35 
obtainable by variation of the parameters 1, 2 is denoted by a point X max. Assum- 
ing the output X is at a point Xo, then movement of the output along a line Zr is 
necessary to attain X max. Zr has components Zl and Z2 in the directions of the 
axes of the parameters 1 and 2 respectively. 

The engine is controlled by a digital computer 28 which is adapted to vary the 
parameters 1 3 2 to maintain X at X max. To this end the computer is programmed 
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to: — 



^ impart perturbation signals Bl, B2 to the respective signals 1, 2. These perturba- 
' tions are added to the signals 1 5 2 at summing junctions 30, 32. 
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(72) 



(71) We 3 Rolls-Royce Limited a British Company of Moor Lane 3 Derby, 
do hereby declare the invention, for which we pray that a patent may be granted to 
us 3 and the method by which it is to be performed, to be particularly described in and 
by the following statement : 

This invention relates to an optimizing control system. 

According to this invention there is provided a system for varying an input of a 
manufacturing or power plant automatically to drive an output of the plant to a 
maximum or minim um value, comprising a digital computer adapted to act on said 
input and respond to said output and programmed to perform iteratively the opera- 
tions of converting a pseudo-random binary sequence into a perturbation signal where- 
in one of the two types of digit of the sequence produces a ramp signal component 
successive ones of which are of opposite slope and the other type of digit produces 
signal components of uniform level connecting the successive ramp components, apply- 
ing the perturbation signal to the plant input, reading the consequential perturbations 
of the plant output, correlating the perturbation signal and the consequential perturba- 
tions to determine the slope of the input/output relationship, generating a correction 
signal which is a function of said slope relationship and indicative of the sense in which 
the input has to be varied to drive the output to said maximum or minirnum value, 
and applying the correction signal to the input. 

An example of a control system according to this invention will now be described 
with reference to the accompanying drawings wherein : — 

Fig. 1 is a system diagram. 

Fig. 2 is a characteristic of certain signals occuring in the system. 

Fig. 3 shows certain aspects of Fig. 1 in greater detail. 

Fig. 4 is a diagram of pulses occuring in certain signals of the system. 

Fig. 1 shows a plant in the form of a gas turbine engine 10 comprising in flow- 
series a compressor 12, a burner 14, a turbine 16 adapted to drive the compressor, a 
jet tube 18 and a variable area exhaust nozzle 20 supported by the jet tube. The 
burner is connected to a fuel supply controlled by a valve 22 which is operable by 
a signal 1. The nozzle 20 is connected to an actuator 24 for varying the area of the 
nozzle. The actuator 24 is operated by a signal 2, The signals l v 2 are hereinafter also 
referred to as "parameters". 

Variation of the parameters 1, 2 causes changes in the output X of the engine as 
measured by a thrust sensor 26. The variation of X is shown by the characteristic 
in Fig. 2 where it is defined by the area 12X. The maximum or optimum thrust 
obtainable by variation of the parameters 1, 2 is denoted by a point X max. Assum- 
ing the output X is at a point Xo, then movement of the output along a line Zr is 
necessary to attain X max. Zr has components Zl and Z2 in the directions of the 
axes of the parameters 1 and 2 respectively. 

The engine is controlled by a digital computer 28 which is adapted to vary the 
parameters 1, 2 to maintain X at X max. To this end the computer is programmed 
to: — 
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a) impart perturbation signals Bl, B2 to the respective signals 1, 2. These perturba- 
tions are added to the signals 1, 2 at summing junctions 30, 32. 
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b) receive the signal X which now contains perturbations occuring in consequence of 
the signals B1, B2. 

c) analyse the signal X with respect to the signals Bl, B2 to determine the effec: of 

these signals on X. 

g d) generate correction signals Zl, Z2 for the varying the respective parameters 1, 2 5 

m the sense causing change of X towards X max. 

To ensure discrimination between the signals Bl, B^^STilbt 
10 g? ^th VEEd versions of its As is ^^S^on^e^ 

eaEtte number of digits of each sequence and also equals the total number of dif- 
■ felSt sW^eTit is possible to have in the group. The sequences Al to An are ^ 

.ffiH^t whkh in turn, is defined as the process of summing the products of the 

auS« men dfe correlation coefficient is +7. This shows that it is possible to , recog- 
25 rizTany oS of the sequences Al to An uniquely from all the other sequence^ . The 25 
Section process will be understood by adding the corresponding digits of the 
r«ofctivrse?uences. This gives a sequence Ax referred to as the collective noise. If 
' rhe P se?uenct AT 11 correlated with any of the sequences Al to A7, say the sequence 
Jl Km t»e found that the correlation products form a sequence Ax idenucal 1 with 
*n Al This ; means that it is possible to have access to any one of the sequences wluch 
30 constuSe ^sequence Ax. It will be noted that, inAe present «cample the correk- 
tTon^oefficient for two identical sequences is +7. This is so ^« ^S.S^^S 
cWn each have one more "1" than "0". If the number of "0"s and I n * 
SSSScTfrnEe equal, as by removing, one "1", then the correlation coefficient for 
35 t^Xticarsequences'is +6. The corresponding correlation 

sTauences would, in the present example, be zero. Thus, if for any reason the PRBN 
hal an "dd '^O" or "1" and thereby does not correlate to zero, die zero correlation 
5n be eTsured by neglecting the odd «0» or "1". The significance of this will 

40 ^Tor^eTurp-r Stahg the correction signals Zl, Z2 the computer is pro- 40 
grammed to perform the following calculation: — 
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(1) Sj " = (n -n)Bjo| 1 / ^ S9n (Bji) X *' ^ " S9 ° (BjF>) * * P 



1=0 



n-1 



(2) Cj 



-Sgn ( Aj (q + \ )) * X (q + 3 ) 



V i = o 
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3 



Sj 

(3) Zj=— 
Q 

wherein: — 

j=the number (1, 2, 3, . . . n) of the parameter. 
Zj=the correction signal for parameter j. 
5 Sj=the slope of parameter j. 5 

Cj=the curvature (i.e. rate of change of slope) of parameter j. 
Aj=the PRBN sequence selected for the control of parameter j. 

Bj=the signal Aj in ramp form, for the purpose of perturbation of parameter j, and 
as varying to either side of a zero mean. 
10 n=the total number of digits in the PRBN sequence; also the total number of para- 10 
meters j which can be controlled by sequences having n digits. 
i=the number 0, 1, 2, . . . n) of any one of the digits. 

X=the output signal of the plant as perturbed by the signal Bj and as varying to 
either side of a zero mean. 

15 P=a value of i such that Bji is positive. 15 
q=a value of i such that Aj (i+i) is negative. 

The above formulae can easily be written in any well-known programming 
language and be read into the computer to constitute a programme therein. The data 

on which the programme has to act is constituted by the signals Aj, Bj and X. 
20 The significance of the various steps of the above formulae and their interaction 20 

with the engine will now be .described in detail. 

The cycle of . the register 33 is initiated by a start signal 34 which determines 

what is referred to as the sampling period of the system. The signal 34 also initiates 

a timing signal 36 (Fig. 4) which determines the time period of the digits i of the 
25 PRBN sequences. The sequences A, A2 are equivalent to square waves and it is ar- 25 

ranged that the time period of each digit i of the sequence extends over two pulses 

of the signal 36. 

The computer programme may be regarded as comprising four sections 38, 40, 
42 and 44 (Fig. 3) concerned respectively with generating the values Bj, Sj, Cj and 

30 Zj. 30 
In section 38 the original Al is converted into the signal Bl which is a ramp 
signal based on a non-return-to-zero response to the signal Al. This means that the 
signal Al is used to generate ramps 46 successive ones of which are of opposite slope 
so that the signal does not have any sharp rise or fall. To make it possible for the 

35 signal Bl to be an identifiable representation of the signal Al only one of the two 35 
types of digit, "1" or "0", are converted into the ramps 46, the other digits remaining 
in the signal Bl as uniform levels 47. Inasmuch as it is a property of the pseudo- 
random binary number that the two types of digit occur respectively an even and an 
odd number of times, the digits chosen for conversion into ramps are those occuring 

40 an even number of times so that the signal Bl. starts and ends at the same level and 40 
no change of level is necessary when the signal is repeated. As mentioned, the signal 
Bl is connected to the summing junction 30 to impart to the signal 1 a perturbation, 
and it will be clear that this perturbation corresponds to the ramps 46, i.e. the signal 
1 is caused to rise and fall in accordance with the ramps 46. The mean level of the 

45 signal Bl is denoted 0 and corresponds to the level which the signal 1 would have in 45 
the absence of the perturbations. The forming of an arithmetical series in the computer 
in order to generate the ramp signal, and the establishment of the mean level of the 
signal, will be well understood in the programming art. 

It will be recalled that the signal Al has one more "1" than "0". This is done 

50 to ensure mat the signal Bl can join the corresponding signal Bl in the next sampling 50 
period without a vertical step. But this odd "1" will be eliminated later in the pro- 
gramme. 

One advantage of using a non-return-to-zero ramp signal for the perturbation 
of the signal 1 is that the engine is capable of responding to such a signal with 
55 reasonably good fidelity. If the signal were a square wave instead of a ramp, then the 55 
engine would have difficulty in responding accurately to the sharp rises and falls of 
such a square wave; also the engine would be subjected to undue wear. 

Similarly to the signal Bl, the signal B2 is a non-return-to-zero ramp signal 
based on the signal A2 and connected to corresponding perturb the signal 2. 
60 The combined effect of the signals Bl, B2 on the engine causes in the output X 60 
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_ , - n v „ chnwn in Fie 4 Fie. 4 also shows the component signal 
Sft&TSSSS & dfSTSS t^sinaYll would have if applied alone, and a 

corresponding component signal X2 which shows the individual effect of the signal 

H7 The sicnals XI, X2 do not exist as such when the parameters 1, 2 are ^^bed 
Srn^neSfHowever they are shown because the object of the correlate as to 5 

^J^^SnSS£STi^ : &S^ the slope SI of the signal 

XI by correlation of the signals Bl, X in accordance with formula (1). The slope 

S2 of the signal X2 is determined similarly. Taking Bl into the formula, as an example 

it will be seen that the expression Sgn (Bli)xXi gives the product of the sign of the 10 

-R1 nt «ch of the seven digits i from 0 to 6. The expression Sgn (Blp)XXp 
Xbtrtcted nfordex » efinate the odd "1". The summation of these products gives 
Ae Safion coemdent. Dividing the summation by n-1 determines the mean 

change of X, (which is the same as the mean change of X) and the further^division 

, t1l , ,-oduius 0 f Bl when i=0 gives the slope SI. As regards the signal X, it will 15 
5 u^demooi how fo form the mefn of the perturbations of the signal X and present 

£dK^ 

^^^^^^^^^^ 20 

of Xi to produce the sequence X. It will be appreciated that the shape of X 
i<! in fact the same as the shape of X but its mean is zero. . „„„„ M , 

It^m als^be seen that the signal X is equivalent to the signal Ax m the general 
example of chelation given earlier herein above, and that the correlation between Bl, 

n X reads on to the component XI of the signal X just as in the general example the 25 
Sau?nce Al read on » the sequence Alx. However, whereas in the V*^*™*** 
the Sidual sequences to bemadded, to form the sequence Ax, each remained un- 
rhLTed i e each had a gain of 1, in the application of the sequences Al, A2 
S^^tSThSa sequences will have gain? other than 1 as is apparent from the 

30 signal X as shown in Fig. 4. The correlation properties, however, of the signal are not 30 

m ^ U T^S^?S , ^S^ in section 42 of the programme by correlation of 

Aj and X. The expression Sgn (Aj(i+i))xX gives the correlation products 

anH rheir summation gives the correlation coefficient. The term Sgn Aj(q+i)xX _ 
35 Slit* SS3 » eliminate an odd "0". Division by (n+l)/2 gives the mean 35 

excursion of the signal X at the midpoints of the digits i, that is at i+J, which i ex- 
SSl bta of curvature. Multiplication by -2/Bjo» gives an approximation 

" ^^^tf^^^'l-oS^the slope is divided bjojjj. 

40 tore to establish the correction signals Z1Z2 which, 

wtioiis 30 32 to drive the parameters 1, 2 to the point X max. lUe tact tnatrne 
sCeTdivided S the curvature means that where the slope is zero, the correction 
Sal will be zero regardless of any inaccuracies in curvature. It can be shown that 
H Sace llx is g based on a quadratic equation then the transition from die 
l;?t Xo to tlie point X max. can be made by a single correlation routine Otherwise 45 

45 S u mo« * chroutines are necessary to attain X max., and it will be understood of 
™urs?£at the computer operates continually to carry out these routines to restrain 
the engine in the X max. condition. 
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Any number of parameters up to n-1 may be used subject however to any pos- 
sible interaction between the parameters in the plant. 

The main advantages of the invention arise, from the fact that PRBN are used 
in combination with a two-level non-square perturbation, signal. A two-level- signal is 

5 one which rises and falls continuously between a highest and a lowest point and does 5 

not halt for a significant length of time at any one intermediate level. The ramp signal 
is an example of this, and such signals are of course easy to produce and process. A 
non-square wave signal is one which, does not have the sharp steps which characterise 
the square wave form but which is in the form, of a ramp 3 a sign wave or any other 

10 form= which is capable of identifying a digit of the PRBN without sharp steps. Such 10 

waves are referred to as non-square waves or as waves having a* finite slope. The 
significance of this requirement is that a square step makes it difficult to determine 

the slope of the signal X when that slope is at or near zero. To illustrate this point 3 

Fig. 4 includes a signal Xla which is the variation- of the signal XI when the para- 
15 meter 1 is close to a maximum point XI max. on the surface 12X. As in- the case 15 

of the signals XI, X2 and X, so also the signal Xla is defined in*, relation, to a zero 
mean to make possible the establishment of the positive and negative values of the 
signal for the purpose of slope correlation. If the signal Bl had been a square wave 

instead of a non-square wave, then, the signal Xla would have consisted of sharp 
20 spikes 48 and this would have made it difficult or impossible to obtain a significant 20 

measurement of the value of the signal at the top of the spike and, as stated, such 
measurement is necessary for forming the mean value which establishes the zero line. 
Thus the invention makes it possible to attain good accuracy at the very point of the 
signal X, i.e. X max., at which the system, is required to maintain the output of tHe 
25 engine. It is emphasized that this accuracy is obtained from a single two level signal 25 
such as a ramp signal and that there is no need to go more complex signals, for 
example signals having three levels to obtain flattening at the top of said spikes, to 
obtain accuracy. 

Further, the present invention makes it possible to obtain accurate measure- 
30 ment of curvature by the use of a simple two-level signal such as a ramp signal. It 30 

will be clear that if the signal Bj had been a square wave, the signals XI, X2, X 
would have substantially vertical sides 54. This would make it impossible to form 

the vahie of X at f+i in formula (2) unless employing a more complicated signal 
pattern, e.g. a three-level signal. 
35 The computer is provided with switch means (not shown) for being connected 35 

or disconnected to the summing junctions 30, 32 and the sensor 26. When the com- 
puter is not in operation the plant is controlled, e.g. manually, by signals 50, 52 to 
the summing junctions 30, 32. 

WHAT WE CLAIM IS: — 

40 1. System for varying an input of a manufacturing or power plant automatically 40 

to drive an output of the plant to a maximum or minimum value, comprising a digital 
computer adapted to act on said input and respond to said output and programmed 
to perform iteratively the operations of converting a psuedo-random binary sequence 
into a perturbation signal wherein one of the two types of digit of the sequence 

45 produces a ramp signal component successive ones of which are of opposite slope 45 
and the other type of digit produces signal components of the successive uniform 
level connecting ramp components, applying the perturbation signal to the plant input, 
reading the consequential perturbations of the plant output, correlating the perturba- 
tion signal and the consequential perturbations to determine the slope of - the input/ 

50 output relationship, generating a correction signal which is a function of said slope 50 
relationship and indicative of the sense in which the input has to he varied to drive 
the output to said maximum or minimum value, and applying the correction signal to 
the input. 

2. System according, to Claim 1 wherein the computer is programmed so that 

55 it is a type of digit occuring an even number of times in the sequence which is con- 55 
verted into said ramp components. 

3. System according to Claim 1 or Claim 2 wherein for the purpose of determin- 
ing said slope the computer is programmed in respect of each digit of the sequence 
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to correlate the sign of *e ~n whichjhe Ration jgg^"^ 
S^SSiSS^Sl." ^io^SUe^ntial outputs 
ata ri^ g ccoSg^r a ny one of tie I^J-£^^££ffi£ 

tion signal as a function of slope and computer is programmed to cor- 
respectively as poame ajd ™* J* = f™?L TmemediaK between the 

"T^S. according to any -erf*. P~^,«^ iSS rfrSL,"^ 

far as the correlation process is c °°«™l d - manufacturing or p0 wer plant automatically 

to ^^urp^S^ A— ^ 38 dC " 

scribed herein with reference to the accompanying drawings. 

For the Applicant 
G. S. COURTNEY 
Chartered Patent Agent 
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